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We solve the Langevin dynamics of d-dimensional ferromagnetic spherical models with interac- 
tions that decay with distance as r~^ d+,J \ The long time dynamics of correlations and responses 
are studied in detail in the different dynamical regimes and the validity of fluctuation-dissipation 
relations (or its violation) are shown. In particular, we show that the fluctuation-dissipation ratio 
■ X(t + t w ,t w ) is assymptotically a fucntion only of the waiting time t w in the aging regime and that 

X — > as t w —* oo. The results are valid in any finite dimension d and for < a < 2 where short 
range behavior is recovered. 

' We also solve the T = Cahn-Hilliard dynamics of this model (conserved order parameter). An 

^ analysis of the multiscaling behavior of the autocorrelation function is presented. 
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I. INTRODUCTION 



Systems with long range interactions are very common in nature. Some important examples are dipolar systems in 
which the interaction decays with distance as 1 /r 3 Jl]-||] ; charged systems with Coulomb interactions oc llr |Q] ; spin 
glasses characterized by RKKY interactions || ; block copolymers |jj7| and models of structural glasses || , to name 
a few. In spite of their ubiquity these kind of systems receive normally much less attention than short range ones, 
perhaps because of its grater analytical complexity. Power law decaying interactions interpolate between the much 
studied, although not realistic, mean field limit and the strictly local nearest neighbour interactions. Although short 
' range interactions are assumed to be dominant in a great variety of systems, this may not be so in many others, spin 
CT) , glasses being a known and controvertial example. 

In pure systems, the dynamics after a quench in temperature from the high temperature phase to the low ordered one 
proceeds by a slow coarsening of domains characteristic of the low energy excitations of the system. This coarsening 
process is characterized by growth laws which show how the typical lenght scales associated with the domains grow 
with time. Typical growing laws are l(t) oc t x / z with z a dynamic exponent characteristic of the universality class 
of the system. Order parameters, correlation and response functions also show typical scaling behaviour. Exponents 
and scaling functions also depend on the dynamics being with conserved or non conserved order parameter and also 
on the nature of the order parameter. From here on we will discuss only results for the vector order parameter 
r T^ • case (n > 1). In the case of long range forces decaying as r -( d + a ) ^ a calculation of the n — > oo limit of the O(n) 
model H gives a growth law l(t) oc t x l a ', for non conserved order parameter. In the conserved case there appears the 
phenomenon of "multiscaling" with two characteristic lenght scales: l\ oc i 1 /( 2 + cr ) and 1% oc (t/ \nt) 1 ^ 2+a \ Here we 
rederive these results working diretly on the spherical model at finite temperature and for general spatial dimension 
d. We will also derive the scaling forms for the structure factor, two times autocorrelations and associated responses. 
In recent years much attention have been devoted to a possible extension of the fluctuation dissipation relations to 
the case of non equilibrium dynamics. It turns out that, in many cases, the fluctuation- dissipation theorem (FDT) 
may be generalized by introducing the concept of time dependent "effective temperatures" |n| . These are difficult to 
obtain both analytically and also to be measured in experiments. The simplicity of the spherical model permits us to 
analyze the character of the violation of the FDT during the coarsening process and to obtain explicitly the effective 
temperature. 
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It is known that the ferromagnetic spherical model with long range interactions has a phase transition in all 
dimensions d (contrary to the short range model which only has transition for d > 2) In d = 1 a ferromagnetic 
phase is present at finite temperature provided that < a < 1 and also in d = 2 provided that < a < 2. Our results 
concerning the ordering dynamics after a quench from the high temperature phase to the low temperature one are 
valid for every a such that < a < s, where s — d for d < 2 and s — 2 for d > 2. The case of a = 2 will be excluded 
as the systems in d = 1 and d = 2 have no phase transitions in this case. For d > 2 and cr > 2 the system has the 
same critical properties of the model with short range interactions. 



II. THE MODEL 



We consider spherical spin models consisting of N continuous spins Si(t) (i = 1, . . . , N) which satisfy for all times 
t the spherical constraint 

N 

i=i 

The Hamiltonian of the system is given by 

n = -J2 J ( r v) s * s i ( 2 ) 

where the sum runs over all distinct pairs («, j) of spins of a d-dimensional hypercubic lattice and is the distance 
vector between sites i and j. The interactions J(jij) decay as a power law of the distance between a pair of spins in 
the following way: 

J( r »j) = J o * J (3) 



where Ylj runs over all sites j ^= i; d is the space dimension, a > 0, Jo > and J(0) = 0. Let us now introduce the 
Fourier transform of the spin variables Si(t) 

*k(i) = -^£>(*)e ik ^ (4) 

j 

with S-k(i) = s kW- Considering periodic boundary conditions the wave vector is k = [k\, . . . , kj) with fcj = 
2Tmi/L,m = 0,±l,±2,...,±(£/2- l),±(L/2) (i = 1, . . . , d) and N = L d . Now defining 

j(^E J W e ' k " ( 5 ) 



the Hamiltonian (0) is diagonalized, 



H = -E J ( k )I^WI 2 ( 6 ) 



with 1 11 



where the sum is over all lattice vectors 1^0. The critical temperature of this model is given by [pr| 
In the long wavelength scaling limit k « 1 Eq.(^) behaves as Jl2] |, 



2 



J(k) m J (l-Cfc c7 + O(fc 2 )) (9) 

where 

C = (10) 



w a fi rf)(T (£) 

2 1 - g 7r d / 2 r(l - (j/2) 
oT(d/2 + cr/2) 



(11) 



v a is the volume of a unit cell in a d-dimensional Bravais lattice £ and fld,a(^) = |1| ^ +CT ^- This approximation 
is valid for < a < 2. At a = 2 logaritmic corrections are present in d = 2. In d dimensions and for a > 2 the 
leading term in the above expansion is 0(k 2 ) and so the critical properties of the model are those corresponding to 
short range interactions. The dynamical properties in this case are also well known (see e.g. Jl3|,|l4j] ) . The analysis we 
present in this work can be extended to the mean filed regime — d < a < 0. However, several works on related models 
p6| p"8[ suggest that the critical properties (both statical and dynamical) in the whole interval are the same as those 
for the case a — —d (mean field). The dynamics for this last case can be almost trivially solved, showing no coarsening 
effects (for instance, the two-times autocorrelation function always decays exponentially with the difference of times). 
Hence, we will not discuss this further in this paper. 

We will concentrate on the analysis of the correlation and response functions. The two-time autocorrelation function 
is defined as 

C ^ t ') = ^E <s r(iK(0> (12) 

r 

= iE^(i,0 (13) 

k 

with Ck(t, t') = (sk(i)s-k(^'))- Note that Ck(t,t) = C(r, £)e~ lk ' r is the dynamic structure factor while C(r, t) = 
(so(t)s r (t)} is the spatial correlation function. The two-time response function is defined as 



6(s r (t)) 



(14) 

h=0 



N ^ ShM' 

r v 

= ^E G k(t,o as) 



N ■ 

k 



with 



where h r (t) is an inhomogeneous external magnetic field and 



G k (M') = t#- (16) 



ik-r 



III. NON CONSERVED ORDER PARAMETER 

In this case we consider a Langevin dynamics for the spins, 

<m zVsM+m, (17) 

where the Lagrange multiplier z(t) enforces the spherical constraint at each instant and £ is a gaussian white noise 
with moments < >= and < £i(*)£j'(*0 >= 2TSij5(t - t') as usual. 
The Langevin equation in Fourier space reads, 

^ = (j(k)-z(i)KW + £kW (is) 
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with < £ k (i) >= and < (*0 >= 2TS k ,-k'S(t - t'). The formal solution of (|T|) is: 

t+\ tn\ JQtft- V *(t')dt' . f J(k)(i-t')- f* z{t")dt" t /,a 

Skit) = s k (0)e w Jo w + / e Jf £ k (i )rft 

Jo 



From ( |l9| ) and setting t > t' we get 



c k (t,t') 



v/s(*)s(f) 

where Ck(0,0) = ^ r C(r, 0)e~ lk r is the initial autocorrelation and 

S(t) = e a /o 

Again using ( |H^ ) we easily obtain 



C k (0,0)e J ( k )(* + *')+2T /* e J ( k )( t+t '- 2t ")S(t")^" 
Jo 



(19) 



(20) 



G k (M0=e J(k)(t -<V|^. 



From the spherical constraint C(t,t) — 1 V£ we obtain a selfconsistent Volterra equation for the function 5(i): 

S(t) = - X] Ck(0, 0)e 2J ( k ) t + 2T e w ( k )(*- t ')H(t')dt' 

1, L 



(21) 



(22) 



In order to solve this equation for S(i) we must specify the initial correlation C(r,0). For a random initial 
configuration with magnetization m we can choose C(r, 0) = <5 ri o + (1 — 5 r fi)mQ, so C k (0, 0) = (1 — to 2 ,) + mo-/V<5 kj0 . 
We will choose too = which is an interesting case for studying the phase ordering dynamics after a quench from the 
disordered to the low temperature ordered phase |[{|. With this choice C k (0,0) = 1. 

The simplest and most relevant case to study is a quench to zero temperature. It has been shown that scaling 
functions and exponents are the same in the whole low temperature phase and the role of temperature flutuations 
only amounts to a renormalization of the amplitudes |l5| ] . We will show this results emerging cleanly from the exact 
solution of the model at temperatures T < T c . 

For T = we obtain from Eq.(E2h in the thermodynamic limit N — > oo: 



J(k)t 



1 



(2tt) 



dke 2J ^. 



Using (|9j) and making a change of variables u = k a t this integral can be explicitely evaluated: 

e 2J t r i-K a t A 



S(t)a 



t d/cr 



-Jo Cm 



du 



For < a < d and t ^S> 1 the asymptotic behavior is 



with 



A = 



(2n) 



S(t) ~ A- 



= 2 Jot 
t d l' J 



(23) 



(24) 



(25) 



i^ e - JoCu du 



where is the volume of a d-dimensional hypersphere of unit radius. Using this result we obtain for the two time 
autocorrelation: 



C k (M') = 



,J {l-Ck"){t+t') 



.4 



(26) 
(27) 
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In particular we note that the dynamic structure factor 



shows the expected scaling form 



with 



Ck (t) = ±-e J ° Ck "Hi (28) 



C k (i) = L d (t)f(kL(t)) (29) 



L{t)=t 1/a , (30) 



thus recovering the n — ► oo limit of the O(n) model P,pp[|. 

It is known that in the asymptotic coarsening regime a ferromagnet ages, that is, two time correlation and response 
functions depend on both times t and t' explicitly through a non trivial scaling form. In coarsening systems the 
dependence is of the form L(t)/L(t') (t,t' ^> 1), i.e., correlations and responses depend on both times through the 
ratio of the corresponding characteristic lenghts. We will analyze the aging dynamics of the present system with the 
inclusion of finite temperature fluctuations. In doing so we will see that two times scaling laws in the aging regime 
are exactly the same, asymptotically, for every temperature T < T c . At finite temperatures the Volterra equation 
( p2| ) can be solved by Laplace transforming S(t): 

POO 

E(s) = / E(t)e- St dt (31) 



Using that Ck(0,0) = 1 the Laplace transform S(s) of s(t) results: 



5( S ) = ^fc^V) (32) 



1 s-2J(k) 

In order to obtain the scaling behavior of S(s) we have to calculate first the function 



In the thermodynamic limit N — > oo 



Introducing (||) and after some algebra one obtains 0,^5) in the scaling regime fccl: 

if d/2 < a < d 

K c -Koc \ elne if cr = d/2 (35) 
if < a < d/2 

where e = s — 2Jq and K c = K(2J$) = f3 c /2 (see Eq.(||)). From now on one would proceed by different routes 
depending on the value of a. Nevertheless it can be shown that the final results for the scaling behaviors of the 
correlations and responses are the same for all < a < 2. Hence, we will only present here the calculations for 
d/2 < a < d. 



With the help of (B2^,(N) and (35) we find 



jk_a( s -2.7 )*-i 

* iS) =l-f + f( S -2J 0) ^ (36) 

with B = j4|T(1 — d/a)\, r(x) being the gamma function. Working in the low temperature limit (3 3> /3 C and for long 
times we can antitransform the above expression obtaining 

A e 2J °* 
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One can see that the only difference with the T = expression ( p5[ ) is the factor (1 — /3 c //3)~ 2 .The time dependence 
and exponents were not affected by finite temperature. Now we can go back and calculate the two time correlation 
(pp|). It is important to note that for the evaluation of Ck(£, t') we need S(t) for all times and not only for long ones. 
A detailed analysis of the behavior of S(t) for short as well as long times shows that its contribution to the long time 
solution of Ck(i, t') will only be finite for times t > t m i c , t m i c being some microscopic timescale. One then obtains 



C k (t,t') 



,-J Ck"(t+t') 



(1 ~ P/Pc) 2 

A 



(tf) 



2Tt' - 
t' 



e 2J Ck°t'u u -i du 



tmic/*' 



(38) 



The behaviour of the autocorrelation depends on the integral which is a function of k a t'. 

Now we will analyze the two different dynamical regimes k a t 3> 1 (fluctuations inside the domains) and k a t <C 1 
(coarsening regime) and show how this reflects in the fluctuation-dissipation relations. 



A. k a i! 2> 1 : bulk fluctuations or quasi-equilibrium behaviour 

In this case the two time autocorrelation ( p8| ) shows the asymptotic behavior 

C k (t,t , )~e- JoCka ( t+ ^ 



A J Q Ck- U' 



2J Ck"t' 



(39) 



and for long times t ^> t' we have that 



C k (M') 



JnCfc CT Vi 



-J Ck°{t-t') 



(40) 



that is, the dynamics becomes stationary with exponential decay of correlations. It is simple to obtain the two time 
response function (Ell) which gives 



G k (M') 



-J Ck°(t~t') 



Comparing Eqs.rt4Q) and (41) we see that the Fluctuation-Dissipation Theorem: 



(41) 



(42) 



is obeyed in this regime. That is, short wavelength fluctuations k 3> 1/L(t) reflect the (local) equilibrium inside the 
domains. We now analyze the more interesting coarsening regime. 



B. kFi! <C 1 : coarsening, non equilibrium behaviour 

The autocorrelation behaves in this regime as 



C k (t,t')~ e 



-J Ck"(t+t') 



-(tt'Y 



2T ft 



+ 0{k a t') 



A — ' (i-|) V*', 
We see that the only effect of temperature on the structure factor is a renormalization of the amplitude: 



c k (t) 



(1 - /3 C / p) pl/a- -2J Ck"t 



A 



(43) 



(44) 



This has the same scaling as the zero temperature limit (|2q). 

In order to analyze the two time scalings in the aging regime we integrate out all modes k in the autocorrelation: 



C(t,t') = 



1 



(2tt)° 



dkC k (t,t') 



(45) 



G 



This gives 



C(t,f) 



(t + t') d/a 

In the aging regime t 3> t' this simplifies to 



(i-/3 c //?) 2 (tt')- + (my* 



(46) 



C(i,0 = 2%-/y/?) 2 (£p (47) 

which shows the aging scaling typical of coarsening systems with L(t) = t~ valid for general vector models. A similar 
computation for the response function gives 

«('.") = (^0)* m 

It can be readly seen that the Fluctuation-Dissipation Theorem ( |42| ) is not obeyed in this regime. Nevertheless it 
can be extended to this non equilibrium regime defining the so called "fluctuation-dissipation ratio" |icj as 

X(t,t ) = ~QQjj~py~ (50) 

dt> 

This function is particularly interesting in the case of complex systems such as spin glasses or structural glasses 
where it can be associated with physical effective temperatures fiof . Its time dependence encodes information on the 
time scales structure of the system p3. In the case of coarsening systems the response is weak and asymptotically 
goes to zero, signalling the weakness of memory effects in these systems p^pl. This can be explicitely seen in this 
model where at long times 

X(t,t') cxit') 1 -^ ^0 (51) 

We see in the above result that X(t, t') is, asymptotically, a function only of £', as is observed in models of structural 
glasses. This implies that, in particular for fixed t' , the fluctuation-dissipation ratio is a constant. The value of the 
exponent 1 — — also shows that the dynamics becomes faster, and memory effects weaker, as the interactions become 
more long ranged. This means that ferromagnetic domain walls move faster as long range interactions become more 
important. 

IV. CONSERVED ORDER PARAMETER 

We now consider the continuous limit of the spherical model, that is, the spin variables Si(t) are replaced by a field 
s(r,i), r now being a d-dimensional continuous position vector, and where the field is subject for all times t to the 
spherical constraint: 

dv[s{v,t)f = V, 

where the integral is carried out over an hypercube of side L with V = L d . We will restric the analysis to zero 
temperature. In the case of conserved order parameter the dynamics is governed by the Cahn-Hilliard equation: 

-^T = W J^J) (52) 

where the Hamiltonian takes now the form 

H' = - JdrJ dr'J (|r - r'|) s(r, t)s(r' , t) + z{t) J dr [s(r, t)} 2 . (53) 
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Transforming Fourier we arrive at 

^ = -k 2 ^ = k*{[J(k)-z(t)}s*} (54) 
where the asymptotic form ([)]) is assumed for J(k) . The formal solution of Eq. (f54|) is now: 

U\ fn\ k 2 J(k)t-k 2 /*' z(t')dt' 

Sk(i) = s k (0)e Jo (55) 

and the structure factor reads: 

C k (t) = Ck (0)e 2k2j ^ 2k2 ^ (56) 

with g(t) = J* z(t')dt' and considering again Ck(0) = 1. Following ]9| we will obtain the scaling form of the structure 
factor by making the reasonable assumption that Ck(t) will present a maximum as a function of k at some k m {t) and 
at long times it will evolve into a Bragg peak. From the spherical constraint (2ir)~ d J dkCk(t) = 1 it follows that 
Ck m (t) scales with k m as 



Maximizing Ck(t) we obtain 



From Eqs.©, © and (§§) we get 



C km (t)^k~ d . (57) 
2 [Jot - g{t)] 



2 + a J Q Ct 



(58) 



ln/c m JoCat 



whose asymptotic solution for long times is 

Id/ 

k n 



(60) 



_J Ca(2 + cr) t 
We can now reconstruct the structure factor obtaining finally: 

C k {t)=[l d {t)f k/k - ] (61) 

in which 

l(t) = (62) 

and the scaling function 

Mx) = (£±^) x 2 _ 2^2+. (63) 

(T (T 

As in the short range case, the phenomenon of "multiscaling" is also present with long range interactions. In our 
case two characteristic lenght scales show up: 

l(t)ott 1/z with z = 2 + a (64) 

and 

kml{t) *{\h) 1/Z (65) 

Interestingly, this second lenght scale grows more slowly than t and will produce a particular scaling form of the 

two time autocorrelations called "sub-aging". This behaviour has been studied in detail recently in for the O(n) 
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model with short range interactions in the n — > oo limit. Here we generalize the result to systems with power law 
interactions. 

Now it is straitforward to write the solution to (pq) at long times: 



s k (t) = s k (0) exp { -k 2+a J Ct + k 2 



2a 



(2 + a) 



J Ct 



2/ a 



lnt 



(66) 



From this result we obtain the two time autocorrelation following the steps done in the non conserved case. The first 
important point to note is that, similar to what happens with short range interactions, the autocorrelation C(t' + t, if) 
relaxes completely in a time r ~ if: 



C(2t',t') oc (t')^ if ~ 1] 



where 



257 



\ 2/a 2 

1 \ 







< 1 



(67) 



(68) 



and we have set JqC = 1. So the dynamics is considerably faster than in the non conserved case. In the regime when 
t <§C t ! we obtain: 



C(t' + T,t') oc (t')~ 



2 \ 2 /" 



exp 



2/tT 



lnf 



2 + a J (2 + cr)er 2 \ f 



(69) 



It is interesting to compare this result with the corresponding one for short range interactions plfl . First of all, if 
<T^2, the correlation decays algebraically with t', as noted above. For fixed t' the relaxation goes as 



C(t'+T,t') OC e-K^) 2 



(70) 



with6=(f) 2/CT (^)^and 



/hTF 



(71) 



This relaxation time is smaller than if and the particular aging dynamics is called "sub-aging" . The same scaling 
is observed in the short range case plj] . 

Finally we want to point out that, although these results are interesting per se as show complex behaviour emerging 
from a closed analytic solution of the model, it is important to note that the presence of multiscaling in dynamics with 
conserved order parameter is limited to the spherical limit of general n-vector models. Simple scaling is recovered 
when corrections to order 1/n are considered 



V. CONCLUSIONS 



Summarizing, we presented the exact solution for the dynamics of the ferromagnetic spherical model with power 
law decaying interactions in arbitrary dimension, after a quench from infinite temperature into the ordered phase. In 
the case of non-conserved order parameter we analyzed the Langevin dynamics of the model, obtaining the exact long 
time scaling form for correlations and responses at finite temperature. In the case of conserved order parameter we 
analyzed a microscopic version of the Cahn-Hilliard dynamics of the model at zero temperature, obtaining the long 
time scaling form of the two time autocorrelation function. These exact results allow us to check at a microscopic 
level several scaling hypotesis about coarsening dynamics in pure systems, which are currently derived at a coarse 
grained level through phenomenologycal Landau-Ginzburg and Cahn-Hilliard equations. In particular, we recovered 
some know results about the scaling behavior of the n — > oo limit of the 0(n) model, by working directly on the 
spherical model. This allowed us to extend those results by obtainig explicit scaling forms for several correlation and 
response functions. 

One new result worth to mention is that of the Fluctuation-Dissipation ratio Eq.(|5l|). This behavior has been 
observed numerically in several non disordered systems PJ23| and it has been conjectured to be characteristic of 
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systems which do not present replica symmetry breaking j24|, i.e. systems with a single pure state. Up to now, the 
present result is one of the few exact solutions available confirming that conjecture. 

A natural and interesting extension of this work is the study of the dynamics of frustrated systems without disorder, 
like models with competing short-range ferromagnetic and long-range antiferromagnetic (e.g., dipolar) interactions, 
in which simulations have shown a rich dynamical behavior |^||. In fact, there is a whole class of different physical 
systems with these characteristic competing interactions, two examples being charged systems with weak coulomb 
interactions B and a proposed model for the behaviour of structural glasses M . 
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